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SUMMARY 

The lift force on a deformable liquid sphere moving in steady, plane 
Poiseuille-Stokes flow and subjected to an external body force is calculated. 

The results are obtained by seeking a solution to Stokes' equations for the 
motion of the liquids inside and outside the slightly perturbed sphere surface, 
as expansions valid for small values of the ratio of the Weber number to the 
Reynolds number. When the ratio of the drop and external fluid viscosities is 
small, the lift exerted on a neutrally buoyant drop is found to be approximately 
one-tenth of the magnitude of the force reported by Wohl and Rubinow (ref. 1) 
acting on the same drop in unbounded Poiseuille flow in a tiibe. The resultant 
trajectory of the drop is calculated and displayed as a function of the exter- 
nal body force. 


INTRODUCTION 

Understanding the dynamics of a single particle and suspensions of parti- 
cles in slow viscous flow is of fundamental importance in many branches of 
science and technology, such as air pollution, raindrop formation, fluidization 
in the chemical process industry, blood flow, the flow of fiber suspensions in 
paper making, et al. (refs. 2, 3). 

The migration of a single spherical particle across the streamlines of a 
nonuniform creeping flow cannot be explained on the basis of Stokes' equations, 
even in the presence of bounding walls; i.e., a sphere experiences no trans- 
verse force at zero Reynolds number. A transverse force does exist theoreti- 
cally if Inertial forces are taken into account (refs. 4, 5, 6). However, 
the situation is different for flexible particles. Experimental observations 

(i^cf. 7) reveal that at low Reynolds numbers, even when a rigid sphere experi- 
ences a negligible transverse force, neutrally buoyant deforming drops (and 


1493 


flexible solid particles) migrate rapidly across streamlines. This suggests 
that the lift force which produces migration arises from the interaction be- 
tween the particle deformation and the surroianding flow field, rather than from 
an inertial effect. 

Chaffey et al. (ref. 8) considered the problem of a deformable liquid sphere 
in Couette-Stokes flow (linear shear). Assuming the drop to be * close* to the 
plane wall bounding the flow, they found that the effect on the deformed drop 
was to produce a force tending to push the drop away from the wall. This force 
has two failings when used alone as a basis for explaining the migration of drops 
in plane Poiseuille flow. First, it neglects the force due to the interaction 
of the parabolic profile with the resultant deformation of the drop. Secondly, 
it cannot be expected to be valid when the drop is not close to the wall. In 
fact, Kamis & Mason (ref. 9) have shown experimentally that the migration rates 
calculated by Chaffey et al. are significantly larger than those which are ob- 
served. The experimental observations were recorded at a considerable distance 
from the walls relative to the particle size. 

In this paper, the hydrodynamic force arising out of the interaction between 
the incident plane parabolic flow and the sphere deformation is calculated. The 
interaction between the drop and the boundary walls is neglected. The analysis 
herein follows that of reference i wherein the case of Poiseuille flow in a 
tube is considered (see also refs. 10, 11, 12). 

FORMULATION 

The drop surface is defined by 

r = 1 + f(0, (})), (1) 

where (r, 6, <j)) are spherical polar coordinates with pole fixed at the center of 
the ^distorted drop and the axis 0 = 0 along the direction of the undisturbed 
velocity U. We denote the velocity of the fluid exterior to this surface by 
V and the pressure by p. All external lengths, velocities and stresses have been 
non-dimensionalized by a , Uq and l-iUpa respectively, where a is the 
radius of the undeformed drop, Uq is a reference velocity to be specified later, 

and y is the viscosity of the fluid outside the drop. Quantities character- 
izing the interior of the drop are distinguished with a prime. It is assumed 
^ -*• ^ 

that {v,p} and {v',p'} satisfy the Stokes' equations and conditions: 
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Av-Vp = 0, V*v=0, v^Uatr*®®, 

Av’ - Vp* = 0, V • v’ = 0, v’ bounded, 

17^^ = = 0, = v'^, = OT^ at r = 1 + f(0,0), 

ex = ea t' + r7^ + r"^ at r = 1 + f(9,0), 
n n 1 / 


(2y 

(3) 

(4) 


where represents the normal velocity component, the tangential velocity 
vector, the tangential stress vector, the magnitude of the normal stress; 
a = y’/y; R^^ and R 2 are the two principal radii of curvature of the drop surface; 
and the dimensionless parameter e = yU^T where T is the constant surface ten- 
sion associated with the interface between the two viscous media. The pressure 
p’ includes a body force K per unit volume which is assumed to act on the drop 
in the posltive-z direction. Equation (4) is Laplace's formula for the equili- 
brium between the normal stress across the surface and the tension and curvature 
of the surface. 

We seek the solution of equations (2) - (4) as expansions valid for small 
values of c which, upon neglecting terms of 0(e), are of the form 


^ -*■ 

V = Vq + e Vi, 


P = Pq + £ Pl» 


v' = v^ + e vj^, p' = e + P^ + e p]^, 

f = e f, . 


(5) 


The solution in general will depend on the parameter a. We shall see (eq. (6)) 
that the term e p^- , which gives the internal pressure of the drop in its 
spherical shape when U = 0, is needed in order to satisfy condition. (4). The 
deformation f is assumed to be 0(e) so that the drop is spherical when U = 0. 

Upon inserting equations (5) into equations (2) and equating coefficients of each 
power of e in each equation, we find that {v. ,p . ,v’ ,p! } for i = 0, 1 satisfy 
the Stokes* equations, Vp^^^ = 0, = U at r = <», and = 0 at r « 

The boundary conditions (3) - (4) must be examined carefully in order to 
determine the contributions at the various orders of e. It can be shown (ref. 1) 
that the normal and tangential components of the velocity and stress vectors can 
be expressed, through a small rotation, in terms of spherical coordinate com- 
ponents. This small rotation is defined with the aid of two Eulerian angles. 
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which are related to the deformation of the drop. A Taylor expansion about 
r “ 1 of the quantities In equations (3) - (4) is then used to transform to an 
equivalent set of botmdary conditions evaluated on the undeformed sphere surface. 
Finally, taking account of equations (5) yields: 

0(e°), 

pIi = 2/a, (6) 

^Or “ ^Or “ ^Oe = ^00’ ^00 “ ^^0’ 

'^Ore “ “'^6r0’ '^Or0 “ “'^Or0* 

0(e). 



9vQr 9f^ 

^Ir + ^1 ~dT “ = 3 ^ Vqq + cosec0 ^ v^^, 

Vi 0 + (3vQg/8r) = [primes], 

Vj ^0 + “ [primes]. 


at r = 1, 
( 8 ) 


9TQr0 df df 

^lr0 ^1 ~ar “90 ^^Orr “ ^000^ ~ cosec0 '^000“ “ [pnmes], 

9Xor0 9f^ 

'^lr0 ^1 3r §0 ^^000 ''' “§0 ("^Orr " "^000^ 

= a [primes] 


^Orr “ “^Orr ~ cosec0 |g (sin0 -|^) + cosec^0 ^-j) fj^ at r = 1. (9) 

90 

Equations (7) involve only the 0(e*^) terms of the expansion, and equation (9) 
gives the deformation f^^ in terms of the quantities determined from the 0(e*^) 
solution. The fvmction f^^ is determined as a particular solution of equation (9) 
subject to the auxiliary conditions that the volume of the drop remains constant, 
and that the centroid of the drop is chosen to coincide with the origin of the 
coordinate system. 


PLANE PARABOLIC FLOW 


We now consider U to be a plane parabolic flow represented by 

U = 8 + 6x + Y x^, (10) 

with respect to a coordinate system fixed at the centroid of the drop, where 
3, 6 and y constants. In order to relate this to plane Poiseuille flow, let 
the drop be located at an orthogonal distance X = b measured from the mid-plane 
between two parallel walls, one stationary, with the positive X-axis pointing 
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In the direction of the other wall moving in the posltive-z direction with velo- 
city U^. Suppose that the drop is moving with dimensionless velocity 
c = dZ/dt , with respect to a coordinate system fixed in the stationary wall,', 
and in the same direction as the Polseullle flow field. Then the velocity of 
the flow with respect to a coordinate system fixed at the centroid of the drop 


where 


U = 1 - (^) 
^0 


(1 + ^) - c 
^0 


X = b + ax. 


Here the reference velocity Uq = G/2u where bQ is one-half the distance 
between the walls bounding the flow, and G 0 is the constant applied pressure 
gradient. By conqiaring equation (10) with equation (11) it follows that 


2^0 ^0 


, Y = -(a/bQ)". 


Since the 0(e ) velocity field satisfies the boundary condition ^ Uk 
as r -*“ <», where k is the \m±t vector in the z direction, it follows from the 
linearity of Stokes’ equations that the O(e^) solution is easily obtained as the 
sum of the known flows past a liquid sphere in uniform stream, linear shear and 
plane quadratic shear. From these solutions the O(e^) hydrodynamic force 
on the drop is calculated to be ^ 


? = 6TTua U 3 + 3 — y] 

^ ori-1 ^ a+1 J 


This is the drag force acting on the vindeformed drop. The value of 3 (or 

equivalently the velocity c of the drop via eq. (12)) is determined so that the 
A 3 

body force a Kk balances Fq, viz. 

^ " 3a+2 Y - k (3^4.2) » 

2 

where k is the non-dimensional body force given by k = K a / h'Uq. 

The deformation f^ is found from equation (9) to be 
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^11 ■*" ^ 12 * 


hi- 


'12 


19 

24 


-11 


a 


.16 


COS0 (cos6) , 


a- 


40 ^ a+1 


10 

[P^CcosS) - ^ cos20 (cos0)]. 


(15) 
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where P 2 » P^ are associated Legendre polynomials- It is seen that the 

deformation is independent of the body force. 

To determine the 0(e) fields {v^,p^, must solve Stokes’ equations 

subject to the interface boundary conditions (8) at r = 1, in addition to having 


v^ vanish at infinity- 


It is found that the drop does not experience any hydro- 

2 

dynamic force owing to incident linear shear (6x) and quadratic shear (yx ) , 
considered independently. However, their combination gives rise to an additional 
0(e) flow needed to satisfy additional interaction terms in the boundary condi- 
tions (8) . It is this interaction field which produces a transverse force on the 
drop. In particular, define the interaction contributions w- and O- ^ by 


Ir 


= -H w 


Ir 


Ir 


Ir’ 


^ ^ .,(1) + .,(2) ^ 

^Ire lr0 lr0 lr0’ 


(16) 


where and v^^^ represent the linear shear and quadratic shear flows, con- 

sidered independently. There are no contributions arising from the uniform flow 
(3) in equations (16) because to 0(e) uniform flow produces no deformation of the 
drop (eqs. (15)). Substituting equations (16) and analogous expressions for 
the other 0(e) velocity and stress components into conditions (8) together with 
the deformations (15), we find that the interaction field must satisfy, at r = 1: 


w 


+ fr (f 


,( 0 ) 


Ir 3r ""1 "Or 


+ hi ^11 “ [primes] 


^^1 ( 0 ) . 
^^00 + 


!£ i 2 ^( 1 ) . ^hl (2) 

90 ^00 30 00 


+ cosecG (■ 


( 2 ) ^ ( 1 ) 


30 


w. 


( 0 ) 


10 9r 


a. ^2 (1) ^ ^ 

^12 ^00 ^11 ^00 


00 ' 90 ''00 

( 2 ) 

(ifl )= [primes]. 


'^10 I 7 ^00^ + ^12 ^^00^^“ [primes]. 


(17a) 
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0 + -i. (f + f + f ^ 

°lre 3r '^^1 ^Or0 ^ ^12 ^Or0 ^ ^11 ^Or0'' 30 


,(0) 

“^000 


3f 


12 


30 




Ihl (t:<2) 

30 '•^Orr 


-•^00^ 


3fii 

■ “3F 


0f, 


0 +i-(f r^2) + . t,(D)_1!ii 

lr0 3r ^11 Or0 ^12 ^Or0-' 30 


'^000 + IT '^000 “ “ [primes] 

( 2 ) 


3f 


12 ^(1) 


000 


+ cosec9 


^!ll (A2) _ (2). ^ 

30 ^^Orr ^000^ 30 


(t: 


30 

( 1 ) 


000 


^(D) 

000^ 


Orr 
a [primes]. 


(17b) 


The quantity in equations (17) denotes the O(e^) radial velocity component 

corresponding to a unifom stream, linear shear flow and quadratic shear flow 
for i = 0, 1, and 2, respectively, and similarly for the other velocity and 
stress components. 

THE LIFT FORCE MD TRAJECTORY 

The 0(e) interaction field can now be found by first inserting the known 
deformations and O(e^) fields into boundary conditions (17). The solid spherical 
harmonics and their coefficients in Lamb's general solution (ref. 13, p. 595) to 
Stokes’ equations for a sphere are then chosen so that conditions (17) are satis- 
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fled. However, tiie perturbation to the 0(e'^) force experienced by the drop after 
deformation depends only upon one solid spherical harmonic of order ~2 in Lamb’s 
general solution (ref. 14). This 0(e) force can be expressed by 


where 


F^ = -ATTueU^a V(r^P_2), 

P _2 = A ^2 ^ ^ P^(cos 0) cos (j) 


(18) 


in the case of a parabolic flow represented by equation (10) 
equation (14) into the computed value of the coefficient A_ 2 » we obtain the 

?i=||TryeUoa6Y ^ 

n=0 


Upon substituting 


(19) 


where 
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10 


11 


12 


-1- f \ 3 / X \ £ , X6\ fr\t^ XX - 

^2 ■*■ 3 ^cH-1^ (a + j^g) (a - ^ » 

. 1 .3 . . 2^. , 1 ? ^3 3119 2 4233 „ j. 321. 

Wl^ ^ ■*■ 19^ ^4^ ^64 ■ 504 ^ 224 56^ » 

» r-i. -v3 /•„ + /1088i7 2 153439 ^ . 28237. 

“ Wr '•“ ir '‘148960 “ " 223440 “ 15960'' * 


(20) 


The direction of the force Fj^ is orthogonal to the direction of the undis- 
turbed flow at infinity. We note that the factor y appearing in equation (19) 
causes the contributions made by the coefficients F^j^ and "5^2 to vanish when the 
quadratic shear portion of the incident parabolic flow is set equal to zero, 
whereas the body force contribution appearing in the coefficient F^^ remains 
finite and proportional to the body force parameter k. In this case repre- 
sents the transverse force exerted on a drop when the incident flow consists of 
a uniform stream plus a linear shear flow. It is observed from the detailed 

calculation which has been performed that the contributions made by the coeffi- 

*)■ 

dents Fjj and F ^2 to Fj^ arise, respectively, from interactions between (a) the 
O(e^) quadratic shear flow and the linear shear deformation fj^j^,and (b) the O(e^) 
linear shear flow and the quadratic shear deformation ^12* It is seen that the 
interaction between O(e^) uniform stream and f ^^2 does not contribute to the 
0(e) force. 

The sign of the force Fj^ can be determined by Inspection of equation (19) 
after inserting the value of b, given by equation (22). The coefficients (20) 
and their sum are plotted against a in figure 1 for a neutrally buoyant drop, 
when k is zero. For small a, the sum Fj^q + F^^^^ + Fj ^2 negative. In this case 
the migration of the drop is always towards the point of zero velocity gradient, 
viz. X = U^bQ/4UQ. It is observed from equations (20) that, in general, the 
direction of migration depends only on the parameters k/y and a and not upon 
radial position (other than which side of the point where dU/ dX = 0 the drop is) . 

We shall now calculate the trajectory of the drop. Its lateral velocity 
d^/dt is determined by equating the Hadamard & Rybczynski drag force for a 
liquid sphere to F^^. Thus, 


jv o aU 


-% F, 

^0 


2 

E 


F = - 19 ( 0 ± 1 ) 

60 n =0 


In' 


( 21 ) 


1500 



(22) 


Integration of equation (21) yields 

U b. 


b = 


w 0 

4U„ 


(bi- 


U b- 
w 0 

4U. 


•) e 


-t/T 






> 




where b^^ is the Initial position of the drop 


The trajectory of the drop is 
obtained by dividing db/dt, given by equation (21) into dz/dt = c, given by 


equation (14) . 
Z -1 


Upon integration we obtain 
2 oH-1 


0 


e(a/bQ)-^F 


(1 + 


3 30H-2 


k - 


a 


^ + — + -^) In 

O ' OTT ' O' 


2 U U 
a . w 

3a+2 ,2 2U- 

bo 0 16Uo 


b. 


U 


w 


4U. 


0 


U 

w 

4U. 


- i c (^)2 - (^)^] + 
^0 *^0 


n ^0 ^0 J ■ 


(23) 


The ratio U^/Uq may be either positive or negative depending on the sign 

of the applied pressure gradient. In undisturbed flow some layers of the fluid 

will move in a direction opposite to when Uq is sufficiently negative. The 

critical pressure gradient occurs when dU/dX = 0 at X = "b^, i.e., when 

U /U^ = -4. Theoretical curves based on equation (23) are shown in figure 2(a) 
w u 

for = 0, Uq > 0 and in figure 2(b) for = ~4, for the parameter values 

a = 0, e = 0.01, a/b^ = 0.1 and b^/b^ = 0.5, and for various values of the body 

force parameter k. These trajectories show that the drop may move inwards or 

outwards. The direction of migration as t for the various values of k can 

be determined by examination of equations (21) and (14). 

Figure 2(a) (wherein U =0) shows that the drop moves in the positive-z 

w 

direction and approaches the axis midway between the walls (b = 0) asymptotically 
for k > -0.011 (to three decimal places). When k = -0.011, the drop travels 
along the line b ® b^^ in the posltive-z direction. For -2.25 < k < -0.011, the 
drop migrates radially outwards, initially in the positive-z direction. When 
k < -2.25, the trajectories are always in the negative-z direction and radially 
outwards. When the above trajectories are compared with those of a drop in 
Polseuille flow in a tube, depicted in Figure 7 of Wohl and Rubinow (ref. 1), 
it is observed that there is negligible difference except in the case |k| « 1. 
This obseirvatlon is to be expected because when a ® 0 and'a/bQ * 0.1, it is 
found that F = 1.337 + 120k (to three decimal places) from equations (21) and 
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(20). The value of F for circular Folseuille flow is 10.075'*" 120k (to three 
decimal places) . The body force parameter k dominates both expressions unless 
|k| « 1. Figure 2(b) (wherein U^/Uq = —4) displays the drop moving in the 
positive-z direction and approaching the stationary wall (b = asymptotically 

for k < -0.011. For -0.011 < k < 6.75, the trajectories are always in the 
positive-z direction and radially outwards. When k > 6.75 the drop migrates 
radially outwards, initially in the negative-z direction. 

SCOPE OF THE ANALYSIS 

The incident parabolic flow is considered to be unbounded in the sense that 
the secondary effect of the wall on the perturbed flow produced by the drop is 
neglected, i.e., interaction between the drop and the boundary walls is 
neglected. Also, the results herein are not applicable unless the effects of 
0(e) and (nonlinear) inertial effects can be neglected. 
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